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Abstract 

Using rigorous constitutive linearization of second variation introduced in |B] we study weak stability of 
homogeneous deformation of the axially compressed circular cylindrical shell, regarded as a 3-dimensional 
hyperelastic body. We show that such deformation becomes weakly unstable at the citical load that coin- 
cides with value of the bifurcation load in von-Karman-Donnel shell theory. We also show that the linear 
bifurcation modes described by the Koiter circle [11] minimize the second variation asymptotically. The 
key ingredients of our analysis are the asymptoticaly sharp estimates of the Korn constant for cylindrical 
shells and Korn-like inequalities on components of the deformation gradient tensor in cylindrical coordi- 
nates. The notion of buckling equivalence introduced in [6j is developed further and becomes central in 
this work. A link between features of this theory and sensitivity of the critical load to imprefections of 
load and shape is conjectured. 

1 Introduction 

Recent years have seen significant progress in rigorous analysis of dimensionally reduced theories of plates 
and shells based on T-convergence [31 EH HI HH H3] . In the framework of these theories one must postulate the 
scaling of energy and the forces apriori, whereby different scaling assumptions lead to different dimensionally 
reduced plate and shell theories. By contrast [S] has no need for such apriori assumptions, while pursuing a less 
general goal of identifying a critical load at which the trivial branch of equilibria becomes weakly unstable. 
This exclusive targeting of the instability without the attempt to capture post-buckling behavior leads to 
significant technical simplifications in a rigorous analysis of the safe load problem for slender structures. 
Present work builds on the ideas of [BJ and applies them to buckling of a circular cylindrical shell under axial 
compression. 

The asymptotics of the critical buckling load predicted by the sign of the 3-dimensional second variation 
agrees with the classical value in the shell theory [HI HH E] . The displacement variations that minimize the 
normalized second variation of the energy are single Fourier modes, whose wave numbers lie on Koiter's circle 
[llj . Using the notion of buckling equivalence [BJ we link the classical variational problem for the buckling 
load obtained from the shell theory and the rigorous analysis of the sign of the 3-dimensional second variation 
of the non-linear elastic energy. 

The problem of buckling of axially compressed cylindrical shells occupies a special place in engineering. 
Cylindrical shells are light weight structures with superior load carrying capacity in axial direction as compared 
to plates of the same thickness. They are ubiquitous in industry. Yet, the classical theoretical value of the 
buckling load is about 4 to 5 times higher than the one observed in experiments [1]. This is understood as a 
manifestation of the sensitivity of the buckling load to imperfections of load and shape [TJ [T7J, [TH1 [5D] . The 
general understanding of the sensitivity of the buckling load to imperfections is via the bifurcation theory 
applied to von-Karman-Donnell equations [TOj HH [7] . The subcritical nature of the bifurcation [9j [7] , where 
the load drops sharply at the critical load is responsible for the observed discrepancy between the theoretical 
and apparent value of the critical load. It is important to note that the formal asymptotics leading to 
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the theoretical value of the critical load gives no indication of the nature of bifurcation and the resultant 
imperfection sensitivity. 

Our approach to buckling of a circular cylindrical shell, along with the rigorous derivation of the classical 
buckling load, may offer an alternative interpretation of the sensitivity of the critical load to imperfections. 
Our analysis makes apparent three different asymptotics of the buckling load, only one of which (the one with 
the largest critical load) is realized in a perfect uniformly axially compressed cylindrical shell. Imprefections 
of load and shape lead to small perturbations in a trivial branch that are shown to be sufficient to change the 
asymptotics of the buckling load. While the rigorous analysis requires more work and is left to future studies, 
we conjecture that the two other buckling loads that are significantly lower that the classical one could lead 
to a more transparent explanation of imperfection sensitivity. 

Recall that a stable configuration y — y(x), x E C M 3 must be a weak local minimizer of the energy 

/ W(Vy)dx- f (y,t(x))dS(x), 
Jn Jan 

where W(F) is the energy density function of the body and t(x) is the vector of dead load tractions. The 
essential non-linearity of buckling comes from the principle of frame indifference W(RF) = W(F) for all 
R G 5*0(3), combined with the assumption of the absence of residual stress Wf(I) = 0. For slender bodies 
these assumptions are fundamental for the computation of the constitutive linearization, which is based on 
the fact that the stresses are small at every point in the body right up to the buckling point, and therefore, 
the material response can be linearized locally. The impossibility of the geometric linearization due to the 
distributed nature of local rotations is the essential feature of slender bodies. 

The constitutively linearized problem can be viewed as the variational formulation of the linear eigenvalue 
problem at the bifurcation point. As such it permits an extra degree of flexibility, as one can replace one 
variational formulation with an asymptotically equivalent one [BJ. This flexibility is used in this paper to 
simplify and eliminate heavy algebraic calculations for the perfect cylinder. 

This paper is organized as follows. In Section [2] we extend the general theory of buckling developed in [BJ, 
so that it applies to more general 3-dimensional bodies, including cylindrical shells. We define an equivalence 
class of functionals characterizing buckling and give criteria for showing that a pair of functionals belongs 
to the same equivalence class. In Section [3] we discuss the asymptotics of the Korn constant for cylindrical 
shells. The technical details of the proof are in Appendix [A] In Section [4] we prove Korn-like inequalities 
where the linear strain is bounded in terms of the specific components if the displacement gradient. In 
Section [5] we compute the compressive part of the constitutively linearized functional whose destabilizing 
action is ultimately responsible for buckling when the load reaches critical. In Section [6] we apply the general 
theory of buckling from Section [2] to axially compressed perfect cylindrical shells and derive the formula for 
the buckling load, as well as a collection of buckling modes parametrized by points on the Koiter's circle 
PI]- In Section [7] we show that the more realistic but also more technically challenging boundary conditions, 
where displacements are prescribed on the top and bottom boundaries of the shell produce exactly the same 
buckling load. This is achieved by exploiting the massive non-uniqueness of the buckling modes for the the 
prescribed average vertical displacement boundary conditions of Section [6j In Section [8] we comment on the 
possibility of the link between the sensitivity to imperfection of the buckling load of a slender structure and 
the presence of "latent" buckling modes with significantly smaller critical loads. 

2 Buckling of slender structures 

Here we revisit the general theory of buckling developed in [BJ. The theory deals with a sequence of pro- 
gressively slender domains Qh parametrized by a dimensionless parameter h. For example, in the case of 
the cylindrical shell, h is the ratio of cylinder wall thickness to the cylinder radius. We consider a load- 
ing program parametrized by the loading parameter A describing the magnitude of the applied tractions 
t(x;h,X) — Xt h (x) + 0(A 2 ), as A —> 0. Here and below O(-) is understood uniformly in a; £ fi/, and 
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h £ [0, ho]. Two fundamental assumptions need to be made in order for the general theory of buckling to be 
applicable^] 

The first fundamental assumption requires the existence of the family of equilibrium deformations y(x; h, A), 
corresponding to the applied loads t(x; h, A) and satisfying the imposed boundary conditions, for any h £ 
[0, ho] and A £ [0, Ao], where ho > and Ao > are some constants. Such a family of equilibria will be called 
a "trivial branch". Neither uniqueness nor its stability is assumed. 

The second fundamental assumption is the absence of "bending modes" in the trivial branch. Here we 
use the term "bending" loosely to indicate any response in which the strain to stress ratio becomes infinitely 
large as h — > even for small applied stress. Formally we assume that the trivial branch lies uniformly close 
to the linearly elastic response: 



sup \\Vy(x;h,\) - I - \X7u h {x)\\ L a 

0<h<h 



< CX 2 , u h {x) = 



dy(x; h, A) 



dX 



(2.1) 



A=0 



when < A < Ao and the constant C is independent of h. 

In [5] we have defined the notion of the near-flip buckling when for any h £ [0, ho] the trivial branch 
becomes unstable for A > X(h), where X(h) —> 0, as h — > 0. This happens because it becomes energetically 
more advantageous to activate bending modes rather than store more compressive stress. 

In hyperelasticity the trivial branch y(x; h, X) is a critical point of the energy functional 



W{Vy)dx 



(t(x;h,X),y)dS. 



(2.2) 



dn h 



In general we restrict y to an affine subspace of W ,oo (0/i;M ) given by 

y£y(x;h,X) + V£, 



(2.3) 



where V£ is a linear subspace of W 1 '°°(flh',^& 3 ) that contains Wo'°°{^h\ K 3 ) and does not depend on the 
loading parameter A. The given function y(x;h,X) £ W 1,ca (ilh; K 3 ) describes the Dirichlet part of the 
boundary conditions, while the traction vector t(x; h, A) describes the Neumann-parlJ^] 

The equilibrium equations and the boundary conditions satisfied by the trivial branch y{x\ h, A) can be 
written explicitly only in the weak form: 



{W F {Vy{x;h,X)),V<j>)dx- 



(t(x;h,X),<p)dS = 



(2.4) 



for every <j> G Vhi where Vh is a closure of V£ in W 1 ' 2 ^^; E 3 ). Differentiating (2.4) in A at A = 0, which is 
allowed due to (2.1), we obtain 



(l Vu h (x), V4>)dx 



(t h (x),cl))dS = 0, 



where L = Wff(I)- 

The energy density W(F) satisfies two fundamental assumptions: 

(PI) Absence of prestress: Wf(I) = 0; 

(P2) Frame indifference: W(FR) = W(F) for every R € 50(3); 
(P3) Positive semidefiniteness property (Lq£,£) > for any £ £ R 3x3 ; 



<f>£V h 



(2.5) 



^ome relaxation of the uniformity in h assumption might be necessary in order to bring our theory to bear on the question 
of sensitivity of the buckling load to imperfections 

2 The use of a general subspace Vjj permits one to describe loadings in which desired linear combinations of the displacement 
and traction components are prescribed on the boundary. 
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(P4) Non-degeneracy (Lo£,£) = if and only if £ T = — £. 

By the properties (P3)-(P4) of Lq, there exists ot\_ > 0, such that 



(L |,0>«L l^ym| 2 , £sym = ^ + £ T ). (2.6) 



2 

The buckling is detected by the second variation of energy 



S 2 £{d>;h,\) = [ (W FF (Vy(x; h, A))V0, S7cf>)dx. 



The second variation is always non negative, when < A < X(h) and can become negative for some choice of 
the admissible variation <fi g Vh, when A > X(h). It was understood in [6] that this failure of weak stability 
is due to the properties (P1)-(P4) of W(F ) and is intimately related to flip instability in soft device. It was 
shown in [6] that in an "almost soft" device the critical load could be captured (under some assumptions) by 
the constitutively linearized second variation 



S 2 £ cl (4>;h,X) 

where 



( {(L o e(</>),e(0))-AK,V0 T V0)}dx, (2.7) 



<r h (x) = -L e(u h (x)), e(u) = ^(Vu + (Vu) T ). (2.8) 

Observe, that cr h is minus the linear stress in the body. This notation is convenient when we are dealing 
exclusively with compressive stresses. The larger the stress, the more compressive the load. Let 



A h = |0 e V h : J {a h , V0 T V0)cte > fjj 



(2.9) 



be the set of all destabilizing variations (see (2.7)). 



Definition 2.1. We say that the loading is weakly compressive if there exists h > so that Ah ^ 0, for 
all h < hg. 



If the stress is not weakly compressive, then both terms in ( 2.7 1 are non-negative for all <fi and the instability 
is clearly impossible. If the loading is weakly compressive then the two terms in (2.7) have opposite signs 
for all variations 6 Ah- The linearized second variation can then become negative for sufficiently large A. 
However, this does not immediately imply negativity of the second variation S 2 £((p;h, A). To examine the 
sign of the second variation we consider the function 

m*(h,X) = inf * 2 £(0; h, A) 



<t>&A h J Qh (a h , Vcf> T X7cf>)dx 

that measures reserve of stability in the trivial branch. Negative values of m* (h, A) signal instability. The 
functional m*(h, A) is based on the representation of the buckling load as a generalized Korn constant in [S]. 

Remark 2.2. The theory of buckling of slender bodies in J^j is based on the analysis of the function 

m(h,\)= mf — — . 

4>&k ||V0|| 2 

We will see that for the axially compressed cylindrical shell 

[ {* h ,Vct> T h V(t>h)dx = o{\\V<t>h\\ 2 ) (2.10) 
Jn h 
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for the minimizer (ph in m*(h,X(h)). Therefore, 

and the link between 5 2 £(4>;h, A) and 6 2 £ c i(4>; h, A) cannot be ascertained. In this paper we will show that 
replacing m(h,X) with m*(h,X) permits to establish the required link. We also observe that the distinction 
between m(h, A) and m* (h, A) is essentially 3- dimensional, since, as we have shown in J^jj, 



Definition 2.3. We say that X(h) is the buckling load if 

X(h) = inf{A > : m*(h, A) < 0}. 

We say that {<ph G Ah ■ h € (0, ho)} is a buckling mode if 

r <5 2 £(<^;MW) 

lim — , „ — = = 1), 

h-yo X(h) J Qh (o- h ,V<pZV(j> h )dx 

where X(h) is the buckling load. 

Definition 2.4. We will call the loading strongly compressive, if X(h) —> 0, as h —> 0. 



According to (2.7), strongly compressive loads imply existence of variations d> G Ah for which the measure 
of compressiveness 



/ {o- h .y4> T yd>)dx 



is much larger than the measure of stability 

& h (d>) = [ (L o e(0),e(0))dx. 
Jn h 

In particular, the body flh has to be slender in the sense that the Korn constant 

K(l «%^W <2 ' u) 

is infinitesimal K(Vh) —> 0, as h —> 0. 

The notion of the buckling mode in Definition |2.3| suggests an extension of the notion of buckling equiva- 
lence to include buckling modes in addition to buckling loads. 

Definition 2.5. Assume J(h,<p) is a variational functional defined on Bh C Ah- We say that the pair 
(Bh, J(h,ip)) characterizes buckling if the following three conditions are satisfied 



(a) Characterization of the buckling load: 

h-^b X(h) 



X(h) 
lim = 1, 



where X(h) is the buckling load and 



X(h)=M J(h,d>). 
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(b) Characterization of the buckling mode: If <ph € Bh is a buckling mode, then 

J{h,<f> h ) 



lim 

h-yO 



X(h) 



(2-12) 



(c) Faithful representation of the buckling mode: If 4>h € Bh satisfies \2.12 ) then it is a buckling mode. 
We remark that by Definition 2.3 of the buckling mode the pair (Ah,3(h, <p)) characterizes buckling, where 

5 2 £(4>;h,X(h)) 



3(h,<j>) = X(h) + 



We envision two ways in which the analysis of buckling can be simplified. One is the simplification of the 
functional 2{h, <p). The other is replacing the space of all admissible functions Ah with a much smaller space 
Bh- For example, we may want to use a specific ansatz, like the Kirchhoff ansatz in buckling of rods and 
plates. According to Definition 2.5 the simplified functional J(h,<p), restricted to the ansatz Bh will capture 
the asymptotics of the buckling load and at least one buckling mode. It is in principle conceivable, that there 
are other buckling modes, not contained in Bh- However, we believe that such a situation is non-generic. 
Even in this non-generic situation our approach would allow to capture all geometrically distinct buckling 
modes, if one can identify the ansatz Bh for each of them. 

It is an elementary observation that the only requirement we need to place on the ansatz Bh is that it 
must contain a buckling mode. 

Lemma 2.6. Suppose the pair (Bh, J(h, «p)) characterizes buckling. Let Ch C Bh be such that it contains a 
buckling mode. Then the pair (Ch, J(h, </>)) also characterizes buckling. 



Proof. Let 



inf J(h,<p). 

4>tc h 



X(h) = inf J(h,<p), \(h) 

<t>&B h 

Then, clearly, X(h) > X(h). By assumption there exists a buckling mode <ph € Ch C Bh- Therefore, 



X(h) J(h,<b h ) 

lim — < nm — — 

>>^o x(h) h-yO \(h) 

since the pair (Bh, J(h, <$>)) characterizes buckling. Hence 

lim ^ = 1, 
\(h) 



(2-13) 



and part (a) of Definition 2.5 is established. 
If 4>h G Ch C Bh is a buckling mode then 



lim 

h->0 



J(h,<t> h ) 
X(h) 



1. 



since the pair (Bh, J(h, <p)) characterizes buckling. Part (b) now follows from (2.131. 
Finally, if <ph € Ch satisfies 



then, 4>h G Bh and by (2.13) we also have 



lim 



lim 

h^Q 



J(h, 4>h) 
X(h) 

J(h, 4> h ) 
X(h) 



1. 



= 1. 



Therefore, <ph is a buckling mode, since, by assumption the pair (Bh, J(h, <fi)) characterizes buckling. The 
Lemma is proved now. □ 



G 



If we replace the second variation S 2 £ (</>; h, A) with the constitutively linearized second variation (2.7) in 
the functional 3(h, <p), we will obtain the functional 

J nh {(T h ,Vcj) 1 V4>)dx <£h(<P) 

which first appeared in [5] where it was shown that the buckling load can be regarded as a generalized Korn 
constant 

X(h) = inf &(h,<p). (2.15) 
4>eA h 

Unfortunately the sufficient conditions for buckling equivalence established in [6 fail to guarantee the validity 
of constitutive linearization in the case of the axially compressed thin-walled cylinder for reasons explained in 
Remark |2.2| Our next theorem proves buckling equivalence and also shows that the constitutively linearized 
functional &(h, <p) captures the buckling mode as well. 



Theorem 2.7 (Asymptotics of the critical load). Suppose that the Korn constant K(Vh) defined by {2.11) 
satisfies 

lim K(V h ) = 0, lim = 0. (2.16) 

h^a fo->o K(Vh) 



Then the pair (Ah, &(h, <p)) characterizes buckling in the sense of Definition 2.5 

Proof. The theorem is proved by means of the basic estimate, which is a simple modification of the estimates 
in [5] used in the derivation of the formula for 5 2 £ c i(<p; h, A): 



Lemma 2.8. Suppose y(x;h,X) satisfies (2.1) andW(F) has the properties (P1)-(P2). Then 

\5 2 £((t>;h,X)-8 2 £ cl (<f>;h,X)\ < C(A||e(0)|| || V0|| + A 2 || V0|| 2 ). (2.17) 

Proof. According to the frame indifference property (P2), W(F) = W(F T F). Differentiating this formula 
twice we obtain 

(W FF (F)t £) = 4(Wcc(C)(F T £), F T i) + 2(W c {C)et), C = F T F. 

We make the following estimate 

I (W CC (C) (F T , F T £) - (Wcc(m €) | < I (Wcc{C) (F T - 1% (F T - | + 

\((W CC (C) - Wcc{T))t,t)\ + 2\(W CC (C)$, (F T - 7)01 
When F is uniformly bounded we obtain 

\(Wcc(C)(F T $,F T t) - (WccWte)] <C(\F- I\ 2 \£\ 2 + \C- I||£ sym | 2 + \F I||£ syni ||£|) . 

Similarly, 

\(w c (c) - w cc (i)(c - /u T oi < c\c - i\ 2 m 2 

When F — Vy(x; h, A) and £ = V</> we obtain, taking into account ( |2.1[ ), that 

\F-I\<CX, \C-I\<CX. 



The estimate (2.17) now follows from the formulas 

4W C c(I) = W FF {I) = L 0) \C — I — 2Xe(u h )\ < CX 2 



□ 
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Let us show that for any e > there exists h(e) > 0, so that for all h < h(e) there exists <fih £ Ah, such 
that 8 2 £(4>h] h, X(h)(l + 2e)) < 0. In that case we will be able to conclude that X(h) < X(h)(l + 2e) for any 
h < h(e). For any fixed e and any h > (for which Ah is non-empty) there exists 4>h £ Ah such that 



& h (<j> h )<X(h)(l + e)£ h (<j> h ), 



(2.18) 



thus, putting A = X e (h) = A(/i)(l + 2e) and 4> — <ph hi (2.17) and utilizing (2.181 we get, 



<5 2 £ (</>„ ; Me (M)< 



+ 



6^(0,0 + C(X(h)\\e(cj>h)\\\\Vcl>h\\ + A(/ l ) 2 ||V0, i || 2 ). 



By (2.6 1 we obtain 



Now, by the Korn inequality we obtain 



e h ((P h ) > a L Je((t> h )\\ 



(2.19) 



6 2 £(cf> h ;h,X £ (h)) < - 



™L ( X(h) X(h) 2 ' 

+ \^KWh) K{V h ) / 



\e{4>h)\\ 



We now see that due to (2.16) we have 5 £{<ph', h, X t {h)) < for sufficiently small h. 

Now, let us show that for any e > there exists h(e) > 0, so that for any h < h(e), any < A < X(h)(l — e) 
and any <f> £ Ah we have 5 2 £(4>; h, A) > 0. Indeed, using (2.17) and the generalized Korn inequality 



©fc(0) > X(h)£ h ((j>) 



we estimate for < A < X(h)(l — e) 



S 2 £(ct>;h, A) > e6 h (<P) - |C|(A(/i)||e(0)||||V0|| + A(/i) 2 || V0|| 2 ) 



Using (2.19) and the Korn inequality we conclude that 



5 2 £{cj);h,X) > ea Lo ~C 



X(h) X(h) 2 
^KWh) + ^yh) 



\e{<j>)\\ 



We now see that 5 2 £(4>;h, A) > for sufficiently small h, which means that X(h) > X(h)(l — e). Therefore, 
part (a) of the theorem is proved. 

We will now establish part (b). Assume now that <ph is a buckling mode. Then ah — > 0, as h — > 0, where 



ah 



5 2 £(4> h -h,X{h)) 
X{h)£ h (<t>h) ' 



Observe that by virtue of (a), the condition (2.16) holds for X(h) replaced by X(h), therefore by (2.16) and 
by the Korn inequality, 



K X(h)\\e(<j> h )\\\\VM < m X(h)\\e(<t> h )\\\\Vct>h\\ < ^ 



X(h) 



i i->o 



&h{4>h) 



and similarly 



ft->o a Lo \\e{<t> h )f 



h-*0 



OLL, 



0. 



(2.20) 



(2.21) 



h^O &h(<j>h) 

Let us substitute = <j>h and A = X{h) into ( 2.17[ ) and divide by &h(4>h)- Using (2.20) and (2.21) we obtain 

X(h) 6 2 £(cf>h;h,X(h)) 



lim 

/i-s-0 



&(h, 4> h ) 



&h{<t>h) 



0. 
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We now use ah to eliminate S 2 £: 



Therefore, 



5 2 £(4>h\ h, X(h)) _ a h X(h) 



& h {4> h ) &(h,4> h ) 

.. (l + a h )X{h) 1 
lim 



fc-s-o R(h, (j> h ) 
Recalling that ah — > 0, as h — > we conclude that 

hm^) 
ft-s-o &(h, 4>h) 

and part (b) follows via part (a). 



Let us prove part (c). Let <ph satisfy (2.12 1. Then, by part (a), 



0, as h — > 0, where 



X(h) 



Therefore, 



\(h)£ h (cl>h) = 



&h{(t>h)- 



Let us substitute <p = 4>h and A = \{h) into ( 2.17[ ) and divide by X(h)€h(<ph) to obtain (using 2.22) 



5 2 £(<p h ;h,X(h)) 



X{h)£ h {4>h) 



fa 



< 



C\l + p h )(X(h)\\e(4>h)\\\\Vc}>h\\+X(h) 2 \\Vcl>h\\ 2 ) 
G h (cj>h) 



(2.22) 



(2.23) 



Note that (2.20) and (2.21) continue to hold for any <ph <E Ah,- Thus (2.23) implies 

S 2 £(<f> h ;h,X(h)) 



lim 

h-to X(h)£ h (<t> h ) 



= 0, 



therefore 4>h is a buckling mode. 



□ 



The passage from the second variation of the trivial branch to the constitutively linearized functional 
R(h, 0) is an important simplification. However, in specific problems we might want to simplify the functional 
fl(h, 4>) even further. For that reason we want to establish general criteria for the validity of replacing one 
functional that characterizes buckling by another. Behind this flexibility is the equivalence relation on the 
space of functionals. 

Definition 2.9. Two pairs (Bh, Ji(h, <p)) and (Bh, Ji{h, 4>)) are called buckling equivalent if the pair 
(Bh, Ji(h, 4>)) characterizes buckling if and only if {Bh, J2Q 1 , <p)) does. 

The notion of B-equivalence was introduced in [S] on the set of functions m(h, A) in order to capture 
buckling load by means of constitutive linearization. Definition |2.9| extends the idea of buckling equivalence 
to functionals in order to capture buckling modes in addition to buckling loads. 

Theorem |2. 10| below gives a convenient criterion of buckling equivalence. 

Theorem 2.10 (Buckling equivalence) . If either 



lim X(h) sup 



1 1 

Ji(h,0) ~ Mh,<f>) 



= 0, 



or 



h^o X(h) <peB h 



(2.24) 
(2.25) 



then the pairs (Bh, J\(h, </>)) and (Bh, J"z(h, 4>)) are buckling equivalent in the sense of Definition 2.9 
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Proof. Let us introduce the notation 



K(h) = inf Ji(h,c/>), i = l,2. 



5\{h) = X(h) sup 



1 



1 



Jx{h,<t>) h{h,4>) 



A W 4>eB h 



Then 



and 



X(h) 




A x (/i) 


A 2 (fc) 




A 2 (ft)| 



sup 



sup 



<*>eB„ Ji{h,<f>) 4>eB H J2{h,4>) 



inf Ji (/i,</»)- inf J 2 (/i,0) 



< W 



A(» A(/i) 

Assume that (Bh, J%(h, (f>)) characterizes buckling. Then we have just proved that if eit her 8 1(h) -4 or 
82(h) — > 0, as h — > 0, then X 2 (h)/X(h) — > 1, as /i — > 0, and condition (a) in the Definition 2.5 is proved for 

J 2 (/*,0). 

Observe that by parts (b) and (c) of Definition |2.5| 0/, € 2?^ is the buckling mode if and only if 



This is equivalent to 



lim Jl j h >* h) = 1. 
AW 



1™ Ji(h, 4>h) 



= 1. 



Therefore, 



since either 



lim = 1, 



h-^o X(h) 
X(h) X(h) 



Ji (h,(j> h ) J 2 (h,cj) h ) 

\Ji(h,<j> h ) - J 2 (h,4>h)\ 
X(h) 



Thus, in view of part (a), <fih is a buckling mode if and only if 

lim J j k ^ h) = 1. 
*-><> x 2 (h) 



□ 



3 Korn's inequality for the perfect cylindrical shell 

Consider the perfect cylindrical shell given in cylindrical coordinates (r, 9, z) as 

C h = I h x T x [0, L], h = [l- h/2, 1 + h/2], 
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where T is a 1-dimensional torus (circle) describing 27r-periodicity in 9. In this paper we consider the axial 
compression of the shell where the displacement <p: Ch — > satisfies one of the following two boundary 
conditions: 

0,(r, 6, 0) = Mr, 0, 0) = Mr, 0, 0) = Mr, 0, L) = Mr, 0, L) = 0, (3.1) 



Mr, 0, 0) = Mr, 0, 0) = Mr, 0, L) = Mr, 0, L) = 0, 



: (r,9,0)d9dr = 0. 



(3.2) 



I h xT 



In the first case the top of the shell is allowed only the vertical displacement and the bottom is kept fixed, 
while in the second case both the top and the bottom of the shell are allowed the vertical displacements. In 
order to eliminate vertical rigid body translations, the average vertical displacement of the bottom edge is 
set to zero. In this paper we will work almost exclusively with the boundary conditions (3.2). In Section [7] 
we will show that our results can be extended to the boundary conditions (3.1). Accordingly, let 

V h = {4>e W 1,a (C h ;R 8 ) : Q holds}. (3.3) 

W h = {(t>£W 1 ' 2 {C h ;R 3 ): Q holds} C V h . (3.4) 
The theorem below establishes the asymptotics of the Korn constant K(Vh)- 
Theorem 3.1. There exists a constant C(L) depending only on L such that 

C(L). 



h^fh 



\e(u)\\ 



(3.5) 



For any u 6 V h . Moreover, K(V h ) = C(L)h 3 / 2 . 

The theorem is proved in Appendix [X] We will also need the following Korn- type inequalities whose proof 
is in Appendix |A"| as well. 



Lemma 3.2. Suppose u r = u$ = at z = and z = L. Then 

||^|| 2 + ll^ 2 || 2 <2|| e N||(||e(w)|| + h r ||), 



(3.6) 



and there exist a constant C{L) > depending only on L and an absolute constant ho > such that for all 
h € (0, ho) and for all L > 



||V«|| 2 <C(L)||eH|| ||e(u)|| + 



(3.7) 



We remark that the power of h in the inequality (3.5) is optimal. Indeed, let nu = [hr 1 '^] (integer part 
of /i -1 / 4 ) and < ?/o < n - Let ip(r),z) be a smooth compactly supported function on (— T]o,r]o) x (0, L). We 
define 

<f) h {6,z) = <p{n h 6,z), 06 [-tt.tt], z€ [0,L]. (3.8) 

Extended 27r-periodically in 9 the function (j> h can be regarded as a smooth function on T x [0,L]. We now 
define the ansatz U h (r, 9, z) as follows 



We compute 



'U?(r,9,z) = -0 h ge (9,z) 
Ltf(r, 9, z) = r$(0, z) + (r - 1)0 *) 
^ U^r, 0, z)= (r- 1)4,00.(9, z) - <f>%(6, z). 



Bm /i 1/4 ||Vl7 h || a = ||^||i, CRa) . 



(3.9) 



(3.10) 
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while 



producing 



limfe 5/4 \\e(U h )f = \\ip,zz\\ 2 L 2 {R2) + ^||¥>,wrolli 2 (R a )> 



K(V h ) = 0(h 3 / 2 ). 
Let (VU) a be the components of Vf/ in cylindrical coordinates 



{a,fi}c{rM,z} 



2 /3- 



We compute 



l(VE7% r | 

II (vu% 



! 1 'Kvc/^ 1 )^!! 2 = o 



= o 



left/ 



h\||2 



|(V(7' i ) I . z || 2 = 0(^ 1 / 4 ) = of^ 



rhM|2 



(v[/%|| 2 + ||(vc/ fc ) ze || 2 = o(/, 3 / 4 ) = o ( l|e( y ) 

||(Vtf h ) M || a + ||(V£/' l U| 2 = 0(/x 5 / 4 ) = O (||e(C/' l )|| 2 ) 



In addition wc also have 



ll^ll 2 = 0{h) = 



|e(L^)| 12 
ft 



(3.11) 



(3.12) 

(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 



Remark 3.3. All functions (f)(r],z) in \3.ty constructed via {3.8) vanish to gethe r with all their derivatives at 
z = 0,L, producing test functions U h that satisfy the boundary conditions (3.1) and thus also, the boundary 
conditions (3.2). Therefore, Theorem 3.1 and Lemma 3.2 hold for the space Wh defined by (3.4). 



4 Korn-like inequalities for gradient components 

In order to understand the buckling of the thin walled cylinders we also need to estimate the L 2 norm of the 
individual components of Vit defined in (3.12) in terms of ||e(it)|| 2 . In this section we will prove that the 
asymptotics ( 3.13 )-( 3.17 ) of gradient components of the test function ( 3.9 1 is optimal. In fact, the inequalities 

2 



(v«MI 



ue.e + u r 



<||e(«)|| 2 , ||(Vw) 



< ||e(«) 



are obvious, while the inequalities 



|(v«M| 2 = 



u r .e - ue 



-7^r l|e(M) " ' 



l(Vu) 



6r\\ — \\ u e,i 



|2<£W|| e(M) ||2 



hVh 



are the immediate consequence of the Korn inequality (3.5). The L norms ||(Vtt) 



and ||(Vit)w 



are within ||e(tt)|| of each other, while the same is true for ||(Vu)g z || = ||ue, z || and ||(Vw) z e||. Thus, in order 
to show that the estimates (3.13)-(3.17) are optimal it suffices to prove the upper bounds on ||ur,z||> Il u e,z|| 
and ||u r ||. 

Lemma 4.1. Suppose u r = u$ = at z — and z = L. Then there exists a constant C(L) depending only on 
L such that 

K|| 2 <^||eNj| 2 . (4.1) 
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K*ll 2 < ^%(u)|| a , 



K,f<^||e(«)f, 



(4.2) 
(4.3) 



Proof. We first observe that the inequality (4.1| follows from ( 4.2 ) via the Poincare inequality, and that (4.3) 
is a direct consequence of (4.1) and (3.6). Hence, we only need to prove (4.2). The proof is based on the 
Fourier series in 9 and z variables. For this purpose we need to extend u(r, 8, z) as a periodic function in 
z£l. Our boundary conditions suggest that u r and ug must be extended as odd 2L-periodic functions, while 
u z will be extended as an even 2L-periodic function. Such an extension results in ff 1 (T 2 x I^M 3 ) functions 
whenever u £ 7? 1 (T x [0, L] x 7^;R 3 ). Here T and T 2 denote 1 and 2-dimensional tori, corresponding to 
27r-periodicity in 8 and [0, 2ir] x [-L, L] -periodicity in (9,z), respectively. Denoting the periodic extension 
without relabeling we have 



u(r,9,z) 



EE' 

m=0 nGZ 



^ n \r,9,z), 



(4.4) 



where 



(/> r (r; m, n) sin 



Uz — <fiz(r; m, n) cos y— — j e , <p z (r;m,n) 
We observe that our periodic extension has the property 



L 

Trmz\ 



e m ", 4>r(r; m, n) = 



1 

ttL 

1 



2-n i-L 



u r sin 



o Jo 

2tt i-L 



irmz 
L 



1 
J 



o Jo 

2tt i-L 



nmz \ 
ug sm ( — - — I e 



L 



j dzd9 



u z cos 



L 





" 1 





~ 




' 1 





" 


Vu(r,6,-z) = - 





1 





Vu(r,9,z) 





1 













-1 










-1 



in cylindrical coordinates, and hence 

i-2tt i-L 



/-2tt f-L 

F tJ \ 2 d9dz = 2 / / , ,, 

J-L JO JO 



FiA 2 d9dz, 



for all cylindrical components Fij of Vm. Therefore, it is sufficient to prove (4.2) for functions of the form 
«<->(r, 9, z) = (Mr) sin (=) , f e (r) sin (=) , f z (r) cos (=)) e** 

Indeed, 

,(™.™))|| 2 = ^||e( w )f. 



m— 1 nG2 



m=0 neZ 



Observe that all functions of the form i/™-™) satisfy the boundary conditions from Theorem 3.1 Therefore, 
Theorem |3 . 1 1 and Lemma|3.2|are applicable to such functions. We now fix m > 1 and n £ Z and for simplicity 



of notation we use (v r , vg,v z ) instead of (vi m ' n \ Vg ,n ' n \ vi" 1 '^). 



We notice that if ||iv|| < 3||e(u)||, then the inequality (3.7) implies that 

IK, 2 || 2 <||v^|| 2 <^|| e w|| 2 , 
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and (4.2) is proved. Let us prove the inequality (4.2) under the assumption that \\v r \\ > 3||e(u)||. In that case 
the inequalities (3.6) and (3.7) become 



IK<?|| 2 + |K*|| 2 < |||e(«)||||«r|| 



||V„|| 2 <^ 



3 

e(«)||||u r 



(4.5) 
(4.6) 



We estimate 



L 2 



n 2 \\v r \\ 2 = \\v r .e\\ 2 < 2\\v r , e - v e \\ 2 + 2\\v e \\ 2 < 5||Vu|| 2 + -—\\ Ve \\* < C(L)\\Vv\\' 



Applying the inequality (4.6) we obtain 



n a IKH<^||e(«)||. 



(4.7) 



We estimate 
and 



Kl| 2 < 2|K + v e ,e\\ 2 + 2|K e || 2 < 5|| e («)|| 2 + 2n 2 \\v e \\ 2 , 

^lhll 2 = K,ll 2 <|l|eWIIIKII, 



due to (4.5). Combining the two inequalities we obtain 



i fi r 2 2 

|K|| 2 <5||e( V )|| 2 + ^^||e(.;)|||K 



(4.8) 



By our assumption ||e(v)|| 2 < ||u r || 2 /9. We use this inequality to estimate the first term on the right-hand 
side of (4.8) and obtain 

Kll < U ^\\e(v)\\. (4.9) 

Using ||e(«)|| < ||u r ||/3 again, we conclude that m < CoL\n\ for some absolute constant Co > 0. In particular, 
n^O, since m > 1. Finally, multiplying now (4.7) and (4.9) we get 

m 2 |M 2 <^%(«)|| 2 , 



which completes the proof. 

Thus, we have established the following Korn-likc inequalities for gradient components. 



□ 



Theorem 4.2. Suppose u r = ug = at z = and z = L. Then there exists a constant C{L) depending only 
on L such that 

||(Vtz) e ^ 2 + ||(V^f <||e(«)|| 2 , (4.10) 



||(V«). re || 2 + ||(V M V|| 2 <^l||e(«)|| 2 , 

IMI 2 + ||(v«U| 2 + IKwui 2 < ^|| e («)|| 2 , 

h 

||(Vu)^|| 2 + ||(Vu) z(5 || 2 < ^P|| e H|| 2 , 



(4.11) 
(4.12) 
(4.13) 



Remark 4.3. Theorem 4-2 is obviously valid for u € Wh C V^. The inequalities (4-l)-(4-3) are also sharp in 
Wh, since, as we mentioned in Remark 3.3 the ansatz (3.9), on which the asymptotically behavior of gradient 
components is achieved contains many functions in Wh- 
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5 Trivial branch in a perfect cylindrical shell 



By perfect cylinder we understand the set, given in cylindrical coordinates as 

C h = {(r,e,z):rel h , 9 e T, z e [0,L]}. 

In order to describe the imposed boundary conditions and loading we need to specify the space Vh and the 
functions t(x;X,h) &ndy(x;\,h) in (2.2) and (2.3), respectively. For the space Vh given by (3.3) we define 



t(x; h, A) 



We also have y(x; A, h) = x + U(x; A, h), where 




e t, z e (0,L), 



o, 

L, 



U v '(r, 9, z- h, A) = a(A)r, U {6 \r, 9, z; h, A) = 



U (z) (r,9,z;h,X) 



0. 



(5.1) 



(5.2) 



where the explicit form of a(A) will be given below for a specific energy satisfying properties (P1)-(P4). 

We observe that during buckling the Cauchy-Green strain tensor C = F T F is close to the identity. 
Therefore, considering the energy which is quadratic in E = (C—I)/2 should capture all the effects associated 
with buckling. Hence, we assume, for the purposes of exhibiting the explicit form of the trivial branch, that 



W(F) 



\(LoE,E), E=\{F T F-I). 



2 v u , „ r 

where the elastic tensor l_o is isotropic. Following Koiter [TT] we consider the trivial branch y(x; h, A) = 
x + U(x) h, A) given in cylindrical coordinates by 

C/ (r) = o(A)r, U {e) = 0, U {z) = -6(A)*, (5.3) 

where the functions a(A) and 6(A) will presently be determined. In cylindrical coordinates we compute 





a 





" 




1 + a 










a+f 








vu = 





a 





, F = 





1 + a 





, E = 


















-6 










1 - 6 










S-6 



Then we compute P = F(LqE), and the traction-free condition Pe r = on the lateral boundary is equivalent 
to the equation 

2a + a 2 = v{2b-b 2 ), 



where v is the Poisson's ratio for l_o. The loading (5.1) implies that (Pe z ,e z ) = —A, which translates in the 
equation 

2!/(2a + a 2 ) + (l-z/)(6 2 - 26) 



£(1-6) 

where E is the Young's modulus. Thus 



2(1 + v)(l - 2v) 
(A) = Vl + K26(A)-6(A) 2 )-1, 



-A, 



where 6(A) is the unique root of Eb{l — b) (2 — 6) = 2A, such that < 6(A) < 1 — l/v3- Such a root exists, 
whenever < A < E/(3\/S). We now see that the fundamental assumption &2.11 is satisfied, since the trivial 
branch parameters do not depend on h explicitly. Choosing A as a loading parameter we obtain 



CTh 





1 



(5.4) 



Remark 5.1. The same trivial branch (5.1), \5.$fy , {5.3) also satisfies the more restrictive fixed bottom 
boundary conditions (3.1). 
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6 Buckling load and buckling mode for the perfect cylindrical shell 



Using the linearized stress (5.4) in the Koiter trivial branch (5.3 1 we compute 

(a h ,Vu T X7u)dx = ||M r ,^|| 2 + \\u z , z \\' 



\ u e,. 



Therefore, the space Ah given by (2.9 1 is simply the set of all functions in Vh, given by (|3.3|) that are not 



independent of z- variable. On the on e han d the estimates (4.2) and (4.3) i mply that R(h, u) > c(L)h, for 
any u G Ah where R(h, u) is given by (2.14). On the other, the test function (3.9) shows that X(h) < C(L)h, 
where A(/i) is given by (2.15). Thus, 

c{L)h<\(h)<C(L)h. (6.1) 

In order to find the exact asymptotics of the buckling load as well as the buckling mode we may simplify 
the functional R(h,u) by observing that ||it r ^|| 2 is much larger than ||m ZjZ || 2 and ||w0, z || 2 , according to the 
estimates (4.2) and (|4.3|, which are asymptotically saturated by (3.9). 



Lemma 6.1. The pair (Ah, &i(h, <j))) characterizes buckling, where 

J Ch (L e(4>),e((f>))dx 



J Ch \4>r,z\ 2 dx 



Proof. By (4.3), (2.19) and (6.1) we have 



1 



1 



L (L o e(0), e(0))da 



< 



C(L) 
Vh 



"to 



C ( 1 

< = o 

- Vh \Hh) 



Therefore by Theorem 2.10 pair (Ah, &i(h, 0)) characterizes buckling. 



□ 



Remark 6.2. Remarks 



boundary conditions (3.1) 



3.3 and 4-3 imply that (6.1) and hence Lemma 6.1 are valid {or the fixed bottom 



6.1 Bounds on the optimal wave numbers 

When l_o is isotropic the minimization of &i(h,d)) can be done in Fourier space. For any function f(r) 
(/ r (r), /g(r), fg(r)) and any m > and n £ Z let 

*m,n(/) = [fr{r) sin [— — j , fg{r) sin [— — j , f z (r) cos y—j—j 

For any m > and n > set 



X(m, n) 



'{Ke(* m ,„(/)):/eC 1 (/ /l ;C 3 )}, n>l 
i* m ,n(f) ■ f e CHh;R 3 ), J Ih h(r)dr = 0, }, n = 0. 



(6.2) 



Observe that X(m, n) c Ah for any integers m > 1 and n > 0, since u € X(m, n) is independent of z if and 
only if to = 0. Let 

\\(h) = inf &i(h,4>), X(h; to, n) = inf &i(h,<j>). (6-3) 



Theorem 6.3. 
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(i) Let \\{h) and X(h;m,n) be given by \ 6.3). Then 

Xi(h) — inf X(h;m,n). 



The infimum in ( 6.4) is attained at m = m(h) and n = n(h) satisfying 

C(L) n(h) 2 C(L) 
m(h) < — — . 



(6.4) 



(6.5) 



x/h ' m(h) ~ y/h 
for some constant C{L) depending only on L. 
(ii) Supposem{h) andn(h) are as in part (i) . Then the pair (X (m(h) , n(h)) , &i(h, </>)) characterizes buckling 



in the sense of Definition 2.5 



Proof. Let us first prove (6.4 1. Let 



a(h) — inf X(h;m, ri). 



1 > 1 

l>0 



ft is clear that X(h;m,n) > Xx(h) for any m > 1 and n > 0, since X(m,n) C Ah- Therefore, a(h) > Ai(/i). 
By definition of a(h) we have 

/ (L e(<t>),e(ct>))dx>a(h)\\cf) r , z f 
Jc h 



(6.6) 



for any <j> G X(m 1 n) and any m > 1 and n > 0. Any </> E -4/i can be expanded in the Fourier series in 6 and z 

oo oo 

4>(r, = ^2^2 ^™-«( r ' 61 ' z )' 

m— n— 

where <pm,n{r, 0, z) € X(m, n). If l_o is isotropic elastic tensor, or even more generally has the form 

(L e, e) = Qi(e rr , e r e,e ge , e zz ) + Q 2 {e r z, e 8z ), 

where Qi(9i, 92, <?3, 94) and (32(91,(72) are arbitrary quadratic forms in their arguments, then the quadratic 
form (Lpe, e) diagonalizes in Fourier space, i.e. 



/ {L e{<j)),e(4>))dx = ^ ^ / (L e(</> m ,„), e(0 m , n ))dx. 

•^ c '» m=0 »=0 

00 00 

ii^ii^EEii^ii 2 - 



We also have 



Therefore, (6.6) implies that 



(Loe(0OT,n) ) e(0 mjn ))(ia; > a 
for every m > 1 and n > 0. Summing up, we obtain that 

(Loe^)^^))^^^^)!!^.,!! 2 

Ch 

for every <p E Ah- It follows that Xi(h) > a(h) and equality is proved. 



m,n) || 2 
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Next we prove (6.5). We observe that, according to Lemma 6.1 



c{L)h < < C(L)h. 

Then \(h;m,n) > Ai(/i) > c(L)h for any m and n. By definition of the infimum, there exist indexes m(h) 
and n(h) such that X(h; m(h), n(h)) < 2C(L)h. By definition of the infimum there exists <p h € X(m(h),n(h)) 
such that &i(h, <f) j < 3C(L)h. Hence, there exists a possibly different constant C(L) (not relabeled), such 
that 

||e(^)f < C(L)/i||# J| 2 = C(L)mW 2 fe||^f . (6.7) 



To prove the first estimate in (6.5) we apply the inequality (3.7| to <p h and estimate ||e(0' l )|| via (6.7). 



Urf = UrJ 2 < WV<t>T < C(L) I m{hfh + 



h\\2 



L 2 



m{h) 
~7h^ 



kh\\2 



Hence 



h- 



m(h)^/h 



>c(L) 



for some constant c(L) > 0, independent of h. Hence, the quantity m(h)vh must stay bounded, as h 
To estimate n{h) we write 



n {hfu h r t = u%t < coiuv^uw 2 + ii^ii 2 ). 
n^ii 2 <Siidi 2 <Sii(v^) 



By the Poincare inequality 
and hence n(/i) 2 ||<^|| 2 < C (L)\\(V <j> h )\\ 2 . Applying and estimating ||e(0' l )|| via MM we obtain 



II 2 

9z\\ , 



n(h) 2 < C(L) (hm(h) 2 + ^jpj , 



from which (6.5)2 follows via (6.5)i. The boundedness of m(h) and n(h) implies that the minimum in (6.4) 
is attained. Part (i) is proved now. 



To prove part (ii) it is sufficient to show, due to Lemma 2.6 that X(m(h), n(h)) contains a buckling mode. 
By definition of the infimum in ( |6.3[ ), for each h e (0, ho) there exists iph € X(m(h),n(h)) C Ah such that 

A x (h) = X(h; m(h), n{h)) < MK^h) <%{h) + (\x{h)) 2 . 

Therefore, 

lim ^ h ^ h) = 1. 
Xi(h) 

Hence, iph G X(m(/i), n(h)) is a buckling mode follows, since the pair (Ah,&i(h,(p)) characterizes buckling. 

□ 



6.2 Linearization in r 

In this section we prove that the buckling load and the buckling mode can be captured by the test functions 
depending linearly on r. In fact we specify an explicit structure that buckling modes should possess. We start 
by defining the "linearization" operator 

C(u) = (v r (0,z),rue(l,0,z) - (r - l)v ri o(0, z), u z (l, 0, z) - (r - l)v r . z (9, z)), 
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where 



v r (9,z) = + u r (r,9,z)dr. 
Jih 



Define the space of vector fields 

X lin = {(f(6, z), rg(6, z) - (r - l)/, fl (0, z),h(6, z) - (r - l)f, z (6, z)) 



(f,g,h)eH 1 (Tx[0,L];R 3 ), f(9, 0) = g(9, 0) = f(9, L) = g(9, L) — 0, / h(9,0)d9 = 0}. (6.8) 

Jo 



2 71 



Incidentally, the test functions (3.9) belong to Xu n . It is also clear that Xu n C Vh- We also observe that if 



u e then C(u) e X{m,n). Let us show that if iph € X (m(h) , n(h)) is a buckling mode then so is 



Theorem 6.4. Suppose m(h) > 1 and > are as m part (%) o/ Theorem 6.3 Let iph € X (m(h) , n(h)) 
be a buckling mode. Then _ 

Mh,^h))<Mh,^h)^+C(L)Vh). (6.9) 

Proof. We will perform linearization in r sequentially. First in u r , then in ug and finally in u z . For this 
purpose we introduce the following operators of "partial linearization" 

mi = C r (u) = (v r (9,z),ue(r,9,z),u z (r,9,z)), 

u 2 = C r .e{u) = (v r (9, z),ru g (l, 9, z) - (r - l)v r> g(9, z),u z (r, 9, z)), 

where 

v r (9, z) = + u r (r,9,z)dr. 
Jin 

For any u € Ah we have 

||e(«i) - e(u)\\ < 2(\\v r _ g - u r j\\ + \\v r>z - u r , z || + \\v r - u r \\). 
By the Poincare inequality we have 

\\v r - w r || 2 < C* /i 2 ||u r ,,.|| 2 < C /i 2 ||e(u)|| 2 . 



Now let u 6 X (m(h) , n{h)) , where m(h) and n(h) satisfy (6.5). Then, 



\\v r ,e - u r ,e\\ = n(h)\\v r - u r \\, \\v r . z - u r . : 



\v r - u r \\. 



(6.10) 
(6.11) 

(6.12) 



Substituting this and (6.11) into (6.10), we get 



||e(«i) - e(u)\\ < C h ( 1 + n(h) + , , [u] 



Taking into account (6.5) we obtain 



|e(tii)-e(u)|| < C(L)Vh\\e(u)\\. 



Therefore 



/ (Loe(ui),e(ui))dB < (l + C(L)Vh) [ {L e(u),e{u))dx. 
Jc h Jc h 



(6.13) 



We now make the next step in the linearization in r and consider u 2 — £r.e{ u )- Observe that e(u 2 ) r g = 0. 
We also see that 

e{u 2 )rr = e(«i)rr, e{u 2 ) zr = e(iti) zr , e(u 2 ) zz = e{u x ) zz . (6-14) 
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The remaining components are estimated as follows 

\\e(u 2 )ee - e(«i)flflll ^ 2 ll4 2 _ u eA> \\e(u 2 )g z - e(«i) fljf || < ||ug - 

Therefore, 



We can estimate 



||e(« 2 )|| 2 < || e ( Ul )|| 2 + Co(||ug - u e , e \\ 2 + - u e , z || 2 ). 
||ug - u^f = n{hf\\u^ - u e f < ^ ||4 2) - n e f, 



due to (6.5 1. Similarly, 



hfi ~ *eA* = '^Pf^\K' M\ z < ^\K> - u e \ 



2 _ ir 2 m(h) 2 (2) 



C(L). 



,(2) 



L 2 



(6.15) 
(6.16) 
(6.17) 

(6.18) 



We now proceed to estimate \\u^ — u g \\. Let 



1 - r 



Therefore, 



;(r, 9, z) = u B .r + v rfi - ug = 2e(ui) r8 H (u r ,e - 



Ml' < 8||e(xM)|| 2 + h 2 \\ Vrfi U8 \\ 2 < 8\\e( Ul )\\ 2 + C(L)Vh\\e( Ul )\\ 2 . 



(2) 

due to the Korn inequality (3.5). Thus, < C(L)||e(ui)||. We can express u e — ug in terms of w as 
follows 

/r nr 
w(t,9, z)dt + / (ug(t, 9,z) — ug(l,9, z))dt. 

Using the Cauchy-Schwartz inequality we have 



Therefore, 



f(t)dt dr<- f(rYdr. 
Ik \Ji J 4 Jlh 



\\ug-uff <-{\\wf + \\ug-u e {i,9,z)f). 



By the Poincare inequality followed by the application of the Korn inequality (3.5) we obtain 

\\ue-ug(l,9,z)\\ 2 <h 2 \\ue >r f<C{L)Vh\\e( Ul )f. 

Therefore, we conclude that 

||^-4 2) l| 2 <^)^l|e(«i)l| 2 - 
\\e(u 2 )\\ 2 <\\e( Ul )\\ 2 (l + C(L)h). 



Hence, (6.16) becomes 



(6.19) 



Recalling (6.14) and (6.15) we get 

||Tr(e( M2 ))-Tr(e(« 1 ))||<C(L)^|| e (« 1 )|| 

Therefore, 



When l_o is isotropic we have 



|Tr(e( M2 ))|| 2 <||Tr(e(« 1 ))|| 2 + C(L)ft|| e (« 1 )r 



(L e(n 2 ),e(« 2 ))d : E = A||Tr(e(M 2 ))|| 2 + 2^||e(« 2 )|| : 



(6.20) 
(6.21) 
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where A and fi are the Lame constants. The inequalities (6.191 and (6.20) imply, using the coercivity of l_o, 

(L e{u 2 ),e(u 2 ))dx<(l + C(L)h) f {L Q e{u x ),e{u x ))dx. (6.22) 

Ch J Cu 

In the last step of linearization we let v = C(u). We compute e(v) rz = and 

e(v) rr = e(u 2 ) rr , e(v) r g = e(u 2 )r8, e(v)gg = e(u 2 )gg. 

We also have 

\\e(v)g z - e(u 2 )g z \\ < 2\\v z .g - u z e\\, \\e(v) zz - e(u 2 ) zz \\ < \\v z . z - u ZlZ \\. 



(6.23) 



Analogously to (6.17) and (6.18) we have 



||e(t>)fc, - e(u 2 )g z f < ^^Wv, - u z f, 



\\e{v) zz - e(u 2 ) zz f < ^^-\\v z - u z f. 



(6.24) 



Integrating the equality u z<r — 2e(u 2 ) rz — v r>z from 1 to r we get 

u z (r, 0,z) - v z = 2 / e(u 2 ) rz (t, 9, z)dt. 



Thus, \\v, - 



< h 2 \\e(u 2 )\\ 2 . Applying this estimate to (6.24) we obtain 



\\e(v) 9z - e(u 2 ) 9z \\ 2 < C(L)h\\e(u 2 )\\ 2 , \\e(v) zz - e{u 2 ) zz \\ 2 < C(L)h\\e(u 2 )\\ 2 
We conclude that 

\\e(v)\\ 2 < \\e(u 2 )\\ 2 (l + C(L)h), ||Tr (e(v))\\ 2 < ||Tr ( e (« 2 ))|| 2 + C(L)fc||e(u a )|| 
and hence for the isotropic and coercive elastic tensor Lq we have 



(L e(v),e(v))dx < (1 + C(L)h) I (L e(« 2 ), e(u 2 ))dx. 

Ch J Ch 



(6.25) 



Finally to prove (6.9) we need to relate kVz an d \\v rz . We estimate «rz|| > \\u r 



J r.z u, r.z\ 



Applying (6.5) and (6.11) to (6.12) we obtain 

K.-«r,.|| <C(L)Vh\\e(u)\\, 



and hence, by (6.11) and (6.5) 



\Vr,z\\ > ll«r,»||-C(L)V^||e(u)||. 



At this point the assumption that u e X (m(h) , n(h)) , where m(h) and n(h) satisfy (6.5) is insufficient 



We also have to assume that u = ip^ e X (m(h) , n(h)) is a buckling mode. Recalling that the pair 
(X(m(h),n(h)),&i(h,cf))) characterizes buckling, we obtain the inequality 



Thus, 



||e(tt)|| 2 < Ofti(h, it) || u r , 2 1| 2 < CX{h)\\u r:Z \\ 2 < Ch\\u r , z \\ 2 . 
hrj > \\u r J(l-C(L)h). 



Combining (6.13), (6.22), (6.25) and (6.26) we obtain (6.9) 



(6.26) 

□ 



Introducing the notation X\i n (m, n) — Xy m n X(m, n) we have the following corollary of Theorem 6.4 
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Corollary 6.5. Let the integers m(h) > 1 and n{h) be as in part (i) of Theorem 6.3 Then the pair 
(Xy ln (m(h), n(h)), &y(h, <f>)) characterizes buckling. 



Proof. By Lemma |2.6| it is sufficient to show that X\\ n (m(h), n{h)) contains a buckling mode. Let tph <= 
X(m(h),n(h)) be a buckling mode. Let us show that C(xph) & Xu n (m(h),n(h)) is also a buckling mode. 



Indeed, by Theorem 6.4 



1 < 



MW< (1 + C (^)Wj), 
Ai(/i) A! (A) 



Taking a limit as h — > and using the fact that iph is a buckling mode, we obtain 

Ai(/i) 

Hence, C(if>h) is also a buckling mode, since the pair (X(m(h),n(h)), Ri(h, cf>)) characterizes buckling. □ 



6.3 Algebraic simplification 

At this point the problem of finding the buckling load and a buckling mode can be stated as follows. We first 
compute 

\\i n (h; m, n) = inf &\(Ji,u). (6.27) 



The we find m(h) and n(h) as minimizers in 

Min(h) = min Xy m (h; to, n). 



(6.28) 



For any u 6 Xy ln (m,n) the integral in r over 1^ can be computed explicitly and the minimization problem 
(6.27) can be reduced to an algebraic problem via the Fourier expansion (4.4). However, the explicit expres- 
sions one obtains are, to use an understatement, unwieldy. Therefore, for the purposes of simplifying the 
algebra, we replace e(u) in the numerator of the functional &i(h, it) by 



E(u) = -(G(u) + G(uf), G(u) 



U r fi - Ug 



U0, r Ug t g + U r Ug iZ 
Uz.r ^z.z 



Let 



$o(h,u) 



J Ch r- 1 (L ^(M),£:(u))da; 



Theorem 6.6. The pair (Ah, &o(h, u)) characterizes buckling. 
Proof. First we observe that 



J Q- - lj (L E(u),E(u))dx < hj (l Q E(u), E(u))dx. 



We easily see that \E(u) — e(u)\ < h\Wu\. Therefore, by the Korn inequality (3.5) we obtain 

\\E(u)-e(u)\\<C(L)h^\\e(u)\\. 

This inequality also implies that 

lleHII^CCL^iie^n + ii^ii, 



(6.29) 
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from which we conclude that ||e(tt)|| < C(L)E(u) for sufficiently small h. Therefore, 

\\E{u)-e{u)\\<C{L)h^\\E{u)l 



(6.30) 



and 



Hence, 



Therefore, 



[(L e(u),e(«)) - (l E(u),E(u))]dx 



< C{L)\\E(u) - e{u)\\\\E{u)\\ < C(L)h^\\E(u)\\ 



f (L e(u),e(u))dx- f r- l (\^E(u),E(u))dx <C{L)h 1 ' A f r- l (L E{u),E{u))da 
Jc h JC h JC h 



It follows that 



X(h) sup 



&o(h,u) &i(h,u) 
1 1 



&o(h, u) &i(h,u) 



< C{L)h 



(6.31) 



i/4 Hh) 



We conclude that condition (2.24) is satisfied, since (Ah,&i(h,u)) characterizes buckling. Then, by Theo- 
rem 



2.10 the pair (Ah, &o(h, it)) characterizes buckling. 



□ 



Remark 6.7. The proof of Theorem 
the fixed bottom boundary conditions 



6.6 



3.1 



uses only the Korn inequality {3.5). Therefore, it is also valid for 



Recalling that X\i n (m(h), n(h)) contains a buckling mode, we have the following corollary of Theorem 6.6 
Corollary 6.8. The pair (Xy ln (m(h),n(h)),&o(h,<p)) characterizes buckling. 

The linearization and passage to the Fourier space make it convenient to introduce the following notation. 

C = {x(r, e,z):r = l, 9 e T, z £ [0, L}} 
is the mid-surface of the undeformed cylindrical shell. For / € W 1,2 (Cq\ M 3 ) we define 



u=U(f), 



For m > 1, n > and f e C 3 we define 



u r = fr(0,z), 

ug = rf e (8, z) - (r - l)/ r ,e(0, z), 
u z = f z (6,z)-(r-l)f r<z (d,z). 



7 r (M)=»* (/,sin(^)e me ) , 



(6.32) 



fg=ne(f e sm(^ 



f z = 5Rc(/ 2 cos(^)e 1 "' 



(6.33) 



We also define 
We compute 



Qo(f) + ^Qx(f) 



B(f) 



B(f) = f \f r . z \ 2 dzd0, 

JCn 
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where fj, is the shear modulus and 

Qo(f) = [ {A\fe,e + /,,, + U\ 2 + 2|/ M + / r | 2 + 2| j^| 2 + |/, ifl + / e ,,| 2 }dx, 

J Co 

Qi(f) = [ {A|/ P ,„ + fr.ee ~ fe,e\ 2 + 2|/r,ee - /e,e| 2 + 2|/ r . zz | 2 + \f e , z - 2f rfiz \ 2 }dx, 

where A = 2v/(l — 2v) and v is the Poisson ratio. The problem of finding the buckling load and buckling 
load is stated as ( 6.27| )-(6.28), where the functional ®.\{h. (ft) is replaced with &o(h, <j>). 

When u = U m _ n {f) G X lin (m, n) the problem (6.27) is purely algebraic and can be solved explicitly. 
However, the minimization in (6.281 is a bit messy. In fact, the functional &o(h, <fi) can be simplified further, 
yielding a very simple algebraic problem for computing the buckling load and the buckling mode. For 
4> = U(f) we define 

Qo(/) + gQ?(/) 



B(f) 



B(f) = I \frA 2 dzd9, 



Q*l(f)= f {Mfr.zz+ frfief + 2\frM 2 +2\fr.zz\ 2 +Mfr, Sz \ 2 }dx, 
JCn 



'Co 

Theorem 6.9. The pair (Xn n (m(h),n(h)),R*(h,c/))) characterizes buckling. 
Proof. We split the proof into a sequence of lemmas. 



Lemma 6.10. Suppose m(h) > 1 and n(h) are integers satisfying (6.5) for all h G (0, ho). Then, There exists 
a constant C{L) > 0, such that for any f G C 3 we have 



Mh,u h )<C(L)&*(h,u h ), 

where u h = U m (K), n (h) (/) € Xu n (m(h) , n(h)) . 

Proof. For simplicity denote ||/|| = ||/||l 2 (c )- If 71 = 0j then 



(6.34) 



A*(h,u) = S$o(h,u) + 



\fr 



from which (6.34) follows. Now, let us assume that n > 1. For each / = J-' mn (f) we have 

QU/) = (A + 2)(m 2 +n 2 ) 2 |/ r | 2 . 

We also have that 

IQiCf) - Qt(f)\ < (A + 2)(2(||/ r) „|| + \\f r ,ee\\)\\he\\ + \\he\\ 2 ) + 411/^,1111/^,11 + ll/^ll 2 , 
Computing the norms in terms of the Fourier coefficients we have 

\Qi(f) - Ql(f)\ < 6(A + 2)(n + l)(m 2 + n 2 )\f e \{\f r \ + |/ 9 |). 

Consider now 2 cases. 

Case 1. \f g \ < 2\f r \. 
In this case we have according to (6.36) and (6.35) that 



(6.35) 



(6.36) 



\Qi(f)-Ql(f)\ 
Qlif) 



< 



36n 
rn 2 + n 2 



< 36, 
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and the inequality (6.34) follows. 

Case 2. \f g \ > 2\f r \. 
Observe that 



thus 



\\fe,8 + fr\\ > \\fe,e\\ IIMI >n\fe\-^> f IM 



Qo(f)>2\\fe,e + M\ 2 > Y \fe\ 2 . 



Dividing (6.36) by this inequality we obtain 



|Qi(/) - Qt(f)\ < 12(A + 2)(m 2 + n 2 )(n+l) 



Thus 



Mf) 



\Sio(h,u)-Sl*(h,u) 



\fr 

\fo 



< 36(A + 2)(m 2 + n 2 ). 



h 2 \Qi(f)-Qt(f)\Qo(f) 



12 



Qo(f) 



B(f) 



< 3/i 2 (A + 2)(m 2 + n 2 )&*(h,u). 



Recalling that m = m(h) and n — n(h) satisfy (|6.5[) we conclude that (6.34) holds 



□ 



Lemma 6.11. Suppose Uh € Xu n (m(h),n(h)) is such that there is a constant C'o independent ofh, such that 
&o(h,Uh) < Coh. Then there exists C± > depending only on L, Lq and Co, such that 



(6.37) 



Proof. Before we start the proof we remark that any buckling mode would satisfy all the conditions of this 
Lemma. Let f h € C 3 be such that u h = U m {h),n{h)(fh)- We also define f h = J" m (ft), n (h)(/fc). We will suppress 
the explicit dependence on h in our notation below, and use m, n, f and / instead of m(h), n(h), fh and fh, 
respectively. 

We start with the application of Lemma |3.2| to Uh € Ah ■ We compute 



\ u e,*\\h(c h ) = h ™ 2 \fe\ 2 + —^-\fe ~ in M 2 > hm 2 \f e \ 2 . 



Then, according to Lemma 3.2 we get 

hm 2 \\f g \\ 2 < 2\\e{u h )\\ L 2 {Ch) (\\u r \\ L 2 {Ch) + \\e{u h )\\ L 2 {Ch) ). 
By coercivity of L and the assumption of the Lemma we have 

C h > ^{Kuh) > - 2 

Ot\_ Q \\ u r,z\\ L 2(c k ) 

Thus, there is a constant C = u\_ Ca such that 

\\<n h )\\ 2 LHCh) < Ch\\ip r , z f L2(Ch) = Ch 2 m 2 \f r \ 2 . 
Using this inequality to eliminate ||e(ztft r )||x,2(c fe ) from the right-hand side of (6.38) we obtain 

\fe\ 2 <c(^ + h\\f r \ 2 <C\f r \ 2 Vh 



(6.38) 



for a possibly different constant C. Using this inequality to eliminate \ fe\ from the right-hand side of (6.36) 
we obtain 

\Qi(f) - QUf)\ < Ch x '\n + l)(rn 2 + n 2 )|/ r | 2 , 
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Recalling the formula (6.35) for Q\{f) we get 

|Qi(/)-Qi(/)l / c^/V + i) 



< 



m 2 + n 2 



It is now clear that 



We also have 



Ch*+iQl(f) 



&*(h,Uh) < Mo(h,u h ) + \&*(h,Uh) - Ro(h,Uh)\. 



Therefore, (6.39) also implies (6.37). 

We are now ready to prove the properties (a)-(c) in Definition |2.5| Let 



(6.39) 



□ 



X*(h)= inf #*(ft,</>)- 

4>eX liD (m(h),n(h)) 



Let iph g X]x n (m(h), n(h)) be a buckling mode, whose existence is guaranteed by the Corollary 6.5 Then by 
Lemma 16.111 we have 

hm — — — = 1. 

h-yO A (ft) 



Part (b) of Definition 2.5 is proved. In particular, we obtain 



«— X*(h) 

lim — ^ < 1. 

h->o A(ft) 



(6.40) 



Now, let 4>h € -^iin(TO(/i), n(/i)) be such that 



lim = 1. 

A* (ft) 



Then by (6.40 ) we have &*(h, 4>h) < Cft for some C > 0, and thus, by Lemma 6.10 we have &o(h, 4>h) < C(L)h. 
The inequality (6.39) then implies that 



Therefore, 



lim ggg^ = 1. 
^0 \*{h) 



Ai(/i) 
hm sup 3v < 1 , 

h^Q A* (ft) 



(6.41) 



which together with (6.40) implies the validity of part (a) of Definition 2.5 In particular, this implies that 

^o(h,4> h ) 
hm — — — — = 1. 
h-to A(ft) 

Hence, <ph must be a buckling mode, since the pair (Xr m (m(/i), n(ft))) characterizes buckling. This proves 
part (c) Definition 2.5 □ 
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6.4 Explicit formulas for buckling load and buckling mode 

In this section we solve the minimization problem 

inf &*(h,u), 

ii£Xij n (m,n) 



for any pair of integers m > 1 and n / satisfying (6.5). Let u £ X\i n (m,n) be given by (6.32). Then 
f r (0,0) = f e (0,O) = f r (e,L) = f g (0,L) = 0, and 



A* (h, u) = fx 



Qo(f) + ^Ql(f) 

B(f) 



B(f) = / \f r ,,\ J dzd0, 

JCo 



where 



Q (f) = f {A\f e ,e + f z , z + f r \ 2 + 2|/ M + f r \ 2 + 2|/,, 2 | 2 + \f M>9 + fe.zfjdx, 

J Co 

01(f) = f {Mfr.zz+ frM 2 + 2\frM 2 +2\fr.zz\ 2 + \Vr,e Z \ 2 }dx. 
JCa 

When f(9, z) is such that u £ Xn m (m,n) we obtain 

Qo($) = h\im4e - fh$ z + </v| 2 + 2\in$ + r | 2 + 2m 2 |0 z | 2 + \in$ x + mfo\ 2 , 
Q{($) = (A + 2)(m 2 + n 2 ) 2 |0 r | 2 . 
The minimum of Qo(</ > ) m (4'8,4'z) is achieved at 

(3A + 4)m 2 + (A + 2)n 2 
(A + 2)(n 2 + m 2 ) 2 ' 
Am 2 - (A + 2)n 2 




(A + 2)(n 2 + m 2 ) 2 ' 
Substituting these values back into the quadratic form Q$ we obtain 

4|0 r | 2 m 4 (A + l) 



and hence 



A* [h; m, n) 



(A + 2)(n 2 +m 2 ) 2 ' 
4m 2 (A+l) h 2 {A + 2)(m 2 + n 2 ) 2 



(A + 2)(n 2 + m 2 ) 2 



Minimizing in (m, n) we obtain 



X* (h) = fi min 



12m 2 



4to 2 (A + 1) h 2 {A + 2)(m 2 + n 2 ) 2 



a" t (A + 2)(n 2 + fh 2 ) 2 
achieved at the Koiter's circle: 



12m 2 



= 2fih 



A+l 



(6.42) 



/i(A + 2)(n 2 + fh 2 ) 2 = 4m 2 v/3(A + 1). 

We see how this equation implies our bounds rh(h) 2 h < C and /m(ft,) 4 < Cm(h) 2 . Hence, for any m 
0,1,.. .,M(h) we define 



(6.43) 
(6.44) 



n(m) 



'2m 



</3(A + l) 
VMA + 2) 



(6.45) 
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Figure 1: Buckling modes corresponding, left to right, to m — 1, m — M(h)/2 and m — M(h) on the Koiter's 
circle. 



where 



M(h) 



2L ^3(A+ 1) 
vr y/h{A + 2) 

The buckling modes can then be labeled by the wave number m = 0, 1, . . . , M(h) and given by 



(6.46) 



4> r = sin(mz) cos(n(m)0), 

, , , , (3A + 4)m 2 + (A + 2)n 2 , r , 

<pg = —hn[m) — = sm(n(mjf/J sin(mz), 



4toV 3 ( A + 1) 

+ 2)n(m) 2 
=^==: cosfr. 

;a + i) 

The figure of the buckling mode corresponding to m = 1 is shown in Figure [T] 



Am 2 - (A + 2)n(m) 2 

/i , cosfmz cos n m )0 ). 

4m v /3(A + l) 



7 Fixed bottom boundary conditions 

If the boundary conditions ( |3.1[ ) are imposed, then we can no longer work with a single Fourier mode space 
X(m,n), since it has a zero intersection with the space Wh defined by (3.4). Hence most of the analysis 
in Section [6] cannot be done for the fixed bottom boundary conditions. However, we can still compute the 
buckling load and exhibit buckling modes by modifying the explicit formulas ( |6.42 1 for the buckling modes 
for the boundary conditions (3.2). According to Remark 6.7 the pair (Ah n Wh,&o) characterizes buckling 
for the boundary conditions ( |3.1| . It is therefore clear that 



Therefore, 



X (h) = inf &o{h, <fi) > Xo(h) = inf Mo(h,cj>). 
4>ew h nA h <PeA h 



lira m > 1. 



h^Q A(/i) 

If we find a specific test function £ Wh H Ah such that 

hm — — — — = 1 
h^o X(h) 
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then 



h^a X(h) 

Which proves that Uh G Wu is a buckling mode and 



1 = Ii m *g£*> > lim >M > 1. 



r a(/i) 



hm — — = 1. 

The idea is to look for the buckling mode in the space XQ n (n) of all functions of the form 



u r = 5Re(0 r (z)e m8 ), cfi r G W 1,2 ([0, L] 

we = 5Re(e m8 (r0 e (z) - (1 - r)inf(z))), <fi e € W 1,2 ([0,i] 
, = Re((^(«) - (1 - r)#.(z))e in9 ), 0, G W^ 2 ([0, L] 



C), #.(0) = 0, 
C), 

C), 0,(0) =0. 



(7.1) 



For any it G -X"{? n (n(/i)) we have 



&o(h, u ) = A* 



12 



(0) 



where 



Qg(0) = / {A|tn& + + <M 2 + 2\ind> 9 + <fi r \ 2 + 2|#| a + |*n0, + ^| 2 }d«, 



Q°(0) = / {A|# - n 2 r - in^fll 2 + 2|nVr + m0 e | 2 + 2|#'| 2 + \<j>' e - 2in<f>' r \ 2 }dz, 
Jo 

B o (0)= / |#| 2 <fe. 
Jo 

Even though the fixed bottom boundary conditions prevent the problem to be diagonalized in the Fourier 
space, it is still useful to represent <f>{z) in the form of Fourier series 



.(z) = E] 4> r (m) sin(mz) , 

m— 1 

oo 

)(z) = 2J ^e(m) sin(mz). 

m— 1 

oo 

;(z) = <fi z (m) cos(mz). 



m=0 



The boundary condition (3.1 1 translate via (7.1) into the additional constraints 



mtfi r (m) = 0, E] <fiz( m ) — 0. 

m—l m— 

In terms of Fourier coefficients we have 

oo oo oo 

QoV) = E = E S O (0) = E 

where 



(7.2) 



m=0 



Q^(0) = A|m0 e - m0 z + </> r | 2 + 2\inxj> e + 4>r\ 2 + 2m 2 \<fi z \ 2 + \in<fi z + m(fi e \ 2 , 
Qln($) = A\m 2 (j) r + n 2 4> r + in$o\ 2 + 2\n 2 % T + in<j> 9 \ 2 + 2m i \<fi r \ 2 + m 2 \<j>e, - 2inc/> r \' 
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<(0) = m 2 |0 r | 2 . 

The fixed bottom boundary conditions do not place any additional constraints on the Fourier coefficients of 
4>q. Therefore, we can minimize Q® m + /i in (j)g to obtain an explicit expression of 4>g{m) in terms of 

cj) r {m) and <j> z (m). However, we may simplify the algebra by recalling that the functional R* could be used 
to compute the buckling load. We therefore determine the relation between 4>g(m) and <j> r (rn), and <f> z (m) by 



minimizing 



We obtain 



tig = m- 



(A + 2)0 r - (A + l)m(j) z 
(A + 2)n 2 + m 2 



(7.3) 



We now cook- up a test function based on (6.42). 
Let m = m(h) be such that 



lim m(h) = oo, lim m(h)V~h = 0. 



(7.4) 



Let n = n(h) = n(m(h)) be given by ( 6.45[ ). We remark that under our assumptions n{h) ^> m{h). Let 

/ sin(mz) sin(m + 2z) 



m + 1 



cos(n9). 



This function satisfies all the required boundary conditions in (|7.1[ ). It's z-derivative also vanishes at the top 
of the cylinder, even though we do not require it. If we define 4> z (m) by (6.42 1 then the resulting function </> 2 
will not vanish exactly at the bottom of the cylindrical shell. That is why we modify (6.42 1 as follows 



where 



<j) z = T(m, n)(cos(m + 2z) — cos(mz)) cos(rt^), 

, (A + 2)n 2 - Am 2 
1 (m, n) — 



(7.5) 



(A + 2)(n 2 + m 2 ) 2 ' 

Once again we observe that <p z vanishes not only at the bottom boundary, but also at the top, accommodating 
even pure displacement boundary conditions on top and bottom edges. We may simplify our test function if 
we retain only the necessary asymptotics as h — > in (7.3) and (7.5): 

sin(mz) sin(m + 2z) 



cos(n6), 

m m + 2 ) 

,-(cos(m + 2z) — cos(mz)) 008(716*), 

n \ _ . 

bg = (7(w, n) sin(mz) — 7(771 + 2, n) sin(m + 2z)) sin(n6 

n 



(7.6) 



Substituting this into ^0 w ^ obtain 



lim 

h->0 



7(777,77) 



,. M /A+l 
lim — ' 
h-tO 2 



Af 



(A + 2)n 2 ■ 

n , (^+2) 2 



m 



(m + 2) 2 



We conclude that 



since 



hm^.lim^2 M 

h->-0 h h-yO h 



A+l 



lim 

h-yO 



(777 + 2) 2 



1. 
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Figure 2: Buckling mode (7.1), (7.6) corresponding, left to right, to m(h) ~ h x / 8 , h x / 4 and h 3 / 8 . 



Thus, the test functions (7.6) are buckling modes for any m{h) satisfying (7.4). Figure [2] shows the 
buckling mode (7.6) for 

'V3(A + 1)V 



m(h) 



h(A + 2) 



a = 1/8, 1/4,3/8. 



8 Discussion 



The key observation in our analysis is that for the test functions U (3.9) we have 



6 h (U h ) = 0(\\e(U h )\\ 2 ) = 0{K{V h )\\VU h \\ 2 ). 
However, the asymptotics of the destabilizing compressiveness term 

€ h {U h )= f (a h ,(VU h ) T VU h )dx 

depends strongly on the structure of the tensor 

er°(0, z) = lim 07, (r, 9, z), 



We saw that for the perfect cylinder cr° = er/,. is given by (5.4) and hence 



£ h (U h ) = O 



\e{U h )\ 2 



If we assume that 



(T h (r, 9, z; h) = cr°(6, z) + hr(9, z) + (r - 1)(t 1 (9, z) + 0{h 2 ). 
Substituting it into the equilibrium equation V • cr h = and passing to the limit as h — > 0, we obtain 



•7) 



rr 




(J 



2 <„ 

n° -A- rr° 

u rz "i - u zz,z 





= 0, 

0. 



= 0, 
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The traction- free boundary conditionon the lateral boundary of the shell r = l±h/2 implies that 

er°{9, z)e r = er 1 ^, z)e r = r(9, z)e r = 



for all (9,z) eTx (0, L). Substituting these equations into (8.8) we obtain 



We also obtain er9 „ = and cr9 



cr rr = U, a r0 = U, a rz = U, <J ee = U. 

e + cr° z 2 = 0. Solving these equations results in the following form for er°: 



<r°(9,z) 





s{9) 
s(9) t{9) - zs'{9) 



(8.9) 



for some functions s{9) and t(9). For generic choices of s(9) and we obtain 

\e(U h )\\ 2 ' 



£h{U' 1 ) = O 



resulting in \{h) — 0(h 5 ^ 4 ). One might conject ure that the imperfections of load can produce a non- 
homogeneous trivial branch leading to <r° given by (8.9 1 and hence to the dramatic change in the asymptotic 
behavior of X(h). 



If we disregard the calculations leading to (8.9) and assume for a moment that <jg g ^ then 

jh\\\2 s 



K{V h ) 

It may be conjectured that imperfection of shape can be mathematically described by such tensor er°. In this 
case the critical load has the asymptotics \(h) ~ K(Vh) = 0(h 3 ^ 2 ). We note that the exponents 5/4 = 1.25 
and 3/2 = 1.5 are close the upper and lower limits of experimentally determined behavior of the buckling 
load 015]. 

Observe that £/i(</>) cannot be larger than \\e(<fi)\\ 2 / K(Vh). Therefore, if the predicted buckling load 
X(h) ~ K(Vh) (Euler buckling in the terminology of 0) then the imperfections of load and shape will have 
negligible effect on the buckling load as in the case of straight solid struts and flat plates. 

Acknowledgments. We are grateful to Lev Truskinovsky for reading the entire manuscript and suggesting 
many improvements in the exposition. This material is based upon work supported by the National Science 
Foundation under Grants No. 1008092. 



A Proof of Theorem 3.1 



The proof of Theorem 3.1 is quite involved and will be split into several relatively simple steps. 



A.l Zero boundary conditionon a rectangle 

For any vector field U = (u, v) on ft — [0, h] x [0, L] and any a E [—1, 1] we define 



G n — 



V x Vy + CM 



e a = -{G a + G T a ) = 



-(Uy+V X ) 



-(Uy + V X ) Vy + 
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Theorem A.l. Suppose that the vector field U = (u,v) G C^f^M 2 ) satisfies u(x,0) = u(x,L) = 0. Then 
for any a G [—1, 1], any h G (0, 1) and any L > 

||Gy 2 < 100||e a || (^ + lkc 

We emphasize that there are no boundary conditions imposed on v(x, y). 
Proof. First we prove several auxiliary lemmas. 

Lemma A. 2. Suppose w(x,y) is harmonic in [0,h] x [0, L], and satisfies w(x,0) — w(x,L) = 0. Then 

Kll 2 < ^IHIIKH + IKH 2 . (A.i) 

Proof. If w(x, y) is harmonic and satisfies w(x, 0) = w(x, L) = then it must have the expansion 



w(x,y) = ^2(A n e™ x + B n e ™ x )sin (^) • 

Therefore, 

2 _ Lh ^ f (imh\ ( 2 ^ d2 „^\ , o,t d \ „/.^ _ sinM^) 



In the expansion of w y we simply multiply A n and i?„ by irn/L, while in the expansion of w x wc multiply 
A„ by irn/L and _B„ by —nn/L: 

m2 Lft, 7r 2 n 2 f / Trnh\ 
IKII =^ 



n— 1 k \ / J 

ll9 L/i 7r 2 n 2 f , / 7rnh\ / , 9 ff nh n9 -Tmh \ n A n 1 



The numbers A„ and £?„ can be arbitrary but such that all the series converge. We can therefore change 
variables 

a n = A n e 2-l , b n = B n e ^ , r„ = 



Then 



h 2 2 ^— ' t 2 L 2 

n=l n 

„n Lh 4^ 7r 2 n 2 



E ^WfrO - l)(«n + bl) + (a n + b n f}, 



|W X || 2 



Obviously, 



n=l 



, E - £2-{(^K) - 1)(« 2 + b 2 n ) + (a„ - M 2 }, 

71=1 



00 2 2 2 2 



IKf- IKf = 2L/iE -£T a nbn < 2LhJ2 ~^-a n b n , 7? = {n G N : a n b n > 0.}. 

n=l n£P 
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Next we estimate 



Similarly, 

llw^H 2 > Lft —^-(4>(t„) - l)a n b n . 



L 2 

riev 



Now we have 



nev ' nev 
Applying the Cauchy-Schwartz inequality we obtain 



E 7 -J]T a n b n < J E —j^T^^n) ~ l) a n&n,JE *( T «) ^2 (^( T «) + i) ™ 6 "' 

neP V nev V ™e"P T ™ 



where 



V>(t) 2 -1 sinh 2 (r)-r 2 ' 
The function $(t) is monotone decreasing on (0, +00), and hence, $(t„) < $(ti) < $(0) = 3. Therefore, 

ll 112 11 112 ^ 2 v /$( 7 rfe/L) 

IPvll - < r IMIIWI, (A. 2) 



and the inequality (A.ll follows. The inequality (A.2) is sharp, since it turns into equality for 

w(x, y) = cosh f~ (x - ~J ^ sin (~ 

□ 

Suppose w(x,y) solves 

(Aw(x,y) = 0, fay) € SI . 
\w(x,y) = u(ar, y), (x,y) € dSl, 

where il = [0, ft] X [0, L], Then Vw is the Helmholtz projection of Vm onto the space of the divergence-free 
fields in L 2 {Sl;R 2 ). 

Lemma A. 3. Suppose that the vector field U = (u,v) £ C 1 (Sl;R 2 ) satisfies U(x,0) = U(x,L) — 0. Let 
w(x,y) be defined by (A. 3). Then for any a € [—1, 1], any h £ (0, 1) and any L > 

llVu- Viy|| < f a/2+ - I ||ej|, Ik - loll < - f V2 + - ) ||e a ||. (A.4) 

\ TT J TT \ TT J 

Proof. This follows your note with tighter constants. 

A(u - w) = Au = (en - 622)2 + 2(ei 2 ) y + csen. 
Multiplying by u — in and integrating we get 

||V(u-w)|| 2 = / {(en ~ e 2 2)(u - w) x + 2e 12 (u - w) y + aen(w - u)}dx 
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By Cauchy-Schwartz inequality we get 

||V(«-«;)|| 2 <||e a ||(V2||V(tt-«;)|| + |a|||tt-HI)- 

By Poincare 

h 2 f h 

\u — w\ 2 dx<^r / \(u — w) y \ 2 dx. 



Hence, 

\\u — w\\ < — ||V(tt — w)\\ 



and (A.4| follows. □ 
We are now ready to prove the theorem. For simplicity of notation we denote 

K = ^(V2+- 



By the triangle inequality and Lemma A. 3 we get 

||G a || 2 = ||e Q || 2 + ^\\v x - Uyf = \\e Q \\ 2 + h\(v x + u y ) - 2{u y - w y ) - 2w y \\ 2 < 

IK|| 2 + 1 1|% + ^|| 2 + 6|K - w y \\ 2 + 6K|| 2 < (4 + 67r 2 K 2 )\\e a \\ 2 + 6\\w y f 



To estimate ||u>y|| via Lemma A. 2 we have by the triangle inequality and Lemma A. 3 

IMI < INI + II"- HI < IMI +K h\\e a \\, 



\w x \\ < \\u x \\ + \\w x - u x \\ < (1 + 7rA )||e Q || 



Therefore, 



Thus, 



where 



a < 2>/3(l + 7rg o; ||u||[[ea[| + (J + + ^ + 



^ 1 1 „ . 1 1 1 1 „ II , /i , \fi , m.n i _M<- Ml- 112 



h 



2 12V3(1 + ttAq) 

GJ < : u e a +Ki\\e c 



Ax = 4 + 6tt 2 A 2 + 6(1 + ttA )(1 + (2^3 + ir)K ). 
Rounding the constants up to the next integer we obtain 

||G Q || 2 <99|| eQ || 2 + -N|||e Q ||. 

The theorem is proved now. □ 

A. 2 Periodic boundary conditions on a rectangle 

Theorem A. 4. Suppose that the vector field U = (u,v) e C 1 ([0, h] x [0, 27r];K 2 ) satisfies u(x, 0) = u(x, 2tt). 
Then there exists an absolute numerical constant Cq > such that for any a G [—1, 1] and any h € (0, 1) 



|G a || 2 <Q,||e Q | 



h 
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Proof. For any fixed t € [0, 2ir] denote Ut — (u — u(a;, t), v + ayu(x, t)) and fit = [0, /i] x [t, t + 2ir]. Observe 
that U t satisfies zero boundary conditions on the horizontal boundary of fl t - We apply now Theorem [A.l| to 
the displacement Ut in fit, 



Note that 



and 



thus 



G a (U t )\\ 2 < ^\\e(G a (U t ))\\ ■ \\u t \\ + 100||e(G a (C/ t ))|| 



G a (U t ) = 



e(G a (U t )) 



h 



v x + ayu x (x,t) Vy + u 



— G a 



-u x (x,t) 
ayu x (x,t) 



u x -u x (x,t) \{u y + v x + ayu x (x,t)) 



\{uy + v x + ayu x (x,t)) 



Vy + U 



= e n 



-u x (x,t) \ayu x {x,t) 



_\ayu x (x,t) 
||G Q || < ||G a (l7 t )|| + (l + 27r)|Ma;,t)||, 



(A.5) 



and 



\e(G a (U t ))\\ < \\e a \\ + (l + 2ir)\\u x (x,t)\\ 



\u t \\ < \\u\\ + \\u(x,t)\\ 



Utilizing now (A.5) and taking into account the last three inequalities we arrive at 

||G Q || 2 <2||G Q ([/ t )|| 2 + 2(l + 27r) 2 ||u :c ( a; ^)|| 2 < 
< ™\\e(G a (U t ))\\ ■ \\u t \\ +2Q0\\e(G a (U t ))\\ 2 +2(1 + 2Tr) 2 \\u x (x,t)\\ 2 < 



< 



200 



z a \\ + (l + 2Tr)\\u x (x,t)\\)(\\u\\ + \\u{x > t)\\) + 



+ 200(||e Q || + (1 + 2Tr)\\u x (x,t)\\) 2 + 2(1 + 2 7 r) 2 |K( : r,t)l| : 



(A.6) 

We complete the proof integrating (A.6 1 in t over [0, L] and then estimating each summand applying the 
Schwartz inequality as follows 

\\u x (x,t)\\dt < (L / \\u x (x,t)\\ 2 dt 2 = 
v Jo 7 

-2 / / 2/ 



L z I I u x (x,t)dxdt) ' = L\\u x \\ < L\\e a \\ 
o Jo 



similarly 



||u(a;,t)||dt<L||u||, 



\\u x (x,t)\\ 2 dt = L\\u x \\ 2 <L\\e a \\ 



u.J.r.l) \ui,: h:\dl •■ ( J \\u x (x,t)\\ 2 dt- J \\u(x, t)\\ 2 dt 



L\\u x \\ \\u\\ < i||e Q |j ||u|| 



□ 
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Let 



U X Uy-V 



V x Vy+U 



e, = -(G„+G?). 



Theorem A. 5. Suppose that the vector field U = (u, v) € C x ([0, h] x [0, 27r];]R 2 ) satisfies U(x, 0) = U(x, 2ir). 
Then there exist absolute numerical constants a > and Co > such that for any h £ (0, a) 



|GU| 2 <Co lie 



Proof. Let V = (u, (1 — x)v), and let 



G 1 = G 1 (V), ei = -{G 1 + Gl). 



We compute 



G, = d 










X 1 


-v ' 




V T 




e 1 = e» + 


X 


2 


V + XV X XVy 




V T 

2 


— XVy 



Thus we immediately obtain that 



|G,|| 2 < ecic^ 2 + IMI 2 + ^(KH 2 + K|| 2 ). 



and 

We also estimate 



Ml < ||G, 



|ei||<l|e,||+/i(IKII + KI|) 



KH < \\vy+u\\ + \\u\\ < ||e* 



Now we apply Theorem A. 4 to the vector field V and a = 1, and obtain 



|G*ir < Co ( ||e 1 || 2 + ||e 1 ||Y + IMr + ^ (Kll + IKII 



Next we apply (A.7| to the terms containing ||ei|| and obtain 



||G*|| < Co ( ||e.||" + H e *ll V + H^KII^II + II^ID + + h dl^ll 



Applying the inequalities (A. 8) to the terms containing ||i> x || and ||tfy|| we obtain 

f ||eJ|M + || u ||||G» 



|GJM<C ||e 



MI 2 + IMI 2 



h 2 \\GJ\ 2 



When h 2 < 1/(2C ) we get the inequality 
l|G,|| 2 <Co(j|eJ 2 

We also have 



MIIIGJI + ||d| 



ColMIHGJI < -IIGJI 2 



r 2 



Thus we obtain 



|G,|| 2 <C ||e 



(A.7) 
(A.8) 



(A.9) 
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To finish the proof of the theorem we write ||m|| 2 using integration by parts and periodic boundary conditions: 



||U|| 2 = (U,U + Vy) + (Uy —V,V) + ||u||' 



Thus, 



u < u e» 



\G. 



Thus, using 2||u||||e*|| < \\u\\ 2 + He*!) 2 we obtain 



\u\\ 2 < ||e s ,|| 2 + 2||G*||||^|| +2||u| 



Applying this inequality to the ||u|| 2 term in (A.9) we obtain 

l|G,|| 2 <C (|| e ,|| 2 + ||e,||M + 

from which the theorem follows. 



IGJIIM 



(A.10) 



□ 



A. 3 Proof of Lemma 3.2 and the Korn inequality 



Step 1. First we prove the analog of Lemma 3.2 in which Vu and e(u) are replaced with 



A = 



U r ,r U r fi - Ug U r ,z 

ue, r ue,e + u r u e ^ z 

U z ,r U z a U z , z 



e(A) = -(A + A<). 



respectively, which is 



\\Ar<C(L)\\e(A)\\ \\e(A)\\ + 



To prove inequality ( A. 11 ) we need to estimate three quantities 

Gi2 = \\u e ,r\\ 2 + \\UrM~Ug\\ 2 1 

In terms of 



G 13 = \\u r 



z,r\\ i 



G 2 3 = \\u z ,g\\ 2 + \\u 8 , z 



E 



i 2 ■ \\u0 ir + Urfi -Ug\\ 



E 13 



E 2 i = G?i 



7^2 r~i2 

^22 — U 



22 - \\u0,e + u r \\ 2 , 



E. 



\\u z , g + U0 g \\ 



f~*2 

^33 — u 33 



(A.11) 



Step 2. In this step we prove the inequality (3.6) in Lemma 3.2 We start with G 23 . Integration by parts, 



using the boundary conditions ug = at z = and z = L and the periodicity in 9 gives 
\(u z ,g,ug ;Z )\ = \(u Z:Z ,ug y g)\ < ||u a , 2 ||||ue,e|| < E 33 (E 22 + IKID- 

Therefore, 

G 2 3 = E 2 23 - 2(«,, fl , u 8jZ ) < E 2 23 + E 2 22 + E 2 3 + 2E 33 \\u r \\ < 2||e(A)||(||e(A)|| + \\u r \\). 



(A.12) 



Step 3. Next we estimate G13. Let us fix 9 e [0, 27r] arbitrarily. Next we apply Theorem A.l| to the 
function 

u(r, z) = {u r (r, 9, z), u z (r, 9, z)) 
and a = 0. We obtain, integrating the inequality over 9 as using the Cauchy-Schwartz for the product term 



G 13 < G E n + E 



2 , E 2 



J X3 
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"'\E U + E 13 + E 33 ) ) < C Q \\e(A)\\ ( \\e(A)\\ + M 



(A.13) 
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where Co is an absolute numerical constant, independent of h and L. 

Step 4. Finally we estimate G\ 2 - Let us fix z e [0, L] arbitrarily. We apply Theorem A. 5 to the function 

u(r, 9) = (u r (r, 9, z), ug(r, 9, z)) 

We obtain, integrating the inequality over z and using the Cauchy-Schwartz for the product term 



Gf 2 <C ||e(A)r + ||e(A)|| 



We estimate via the ID Poincare inequality 



t2 t2 or2 

Ml 2 < -\\uoA 2 < ~G 2 23 < — (HA)\\ 2 + He(A)HIKII). 

7T Z 7T Z 7T Z 



(A.14) 



Thus, there exists a constant C(L) < Co(L 2 (<r + 1) + 1) such that 

G\ 2 <C{L)\\e{A)\\ f||e(A)|| + Jl 



Next we prove the analog of the Korn inequality in which again Vu and e(u) are replaced with A and e(A) 
respectively. Integrating the inequality (A. 10) in z and using Cauchy-Schwartz for the product term we obtain 



Ml 2 < l|e(A)|| 2 + 2|| A|||| Ufl || + 2\\u e \\ 2 < ||e(A)|| 2 + e 2 ||A|| 



1 



U0\ 



for any e > 0. The small parameter e £ (0, 1) will be chosen later in an asymptotically optimal way. Applying 

L 2 



(A.14) we obtain for sufficiently small e 



Therefore, 



Thus, 



Kf < ( ~ + 1 ) l| e (A)|| 2 + 6 2 ||A|| 2 + ^-||e(A)|||K||. 



L 2 



\uJ\ < 2 



< V2 



1 || e (A)||^ + 2 £ 2 ||A|| 



+ 1 ||e(A)||+e||A|| 



Substituting this inequality to (A.ll) we conclude that there is a constant C(L), depending only on L such 
that 

\\A\\ 2 <C(L)(-^ + ^) \\e(A)\\ 2 . 



he 2 h 2 



We now choose e = h 1 ^ to minimize the bound: 



| A ||2 < C ( L )^, ^ .11-! 



hy/h 



HA) 



Theorem 3.1 is now an immediate consequence of (A. 15) and the obvious inequality 

\\e(U) - e(A)\\ 2 < ||VC/- A|| 2 < h 2 \\A\\ 2 . 



Observe that we get from (A. 16) and (A. 15) that, 

\\e(U) - e(A)|| < h\\A\\ < C(L)/i 1/4 ||e(A) 
Lemma 3.2 will now follow from ( A.16[ ), (A. 17) and (A.ll). 



(A.15) 

(A.16) 
(A.17) 
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